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Abstract 

■ A few simply-stated rules govern the entanglement patterns that can occur in mutually unbiased 

basis sets (MUBs), and constrain the combinations of such patterns that can coexist (ie, the 

Q_il stoichiometry) in full complements of {p + 1) MUBs. We consider Hilbert spaces of prime power 

' dimension (as realized by systems of N prime-state particles, or qupits), where full complements 

^ ' are known to exist, and we assume only that MUBs are eigenbases of generalized Pauli operators, 

without using a particular construction. The general rules include the following: 1) In any MUB, 

^ I a particular qupit appears either in a pure state, or totally entangled, and 2) in any full MUB 

OO ■ 

complement, each qupit is pure in (p+l) bases (not necessarily the same ones), and totally entangled 

00 ; 

. in the remaining (p — p). It follows that the maximum number of product bases is p + 1, and 

\ when this number is realized, all remaining {p^ — p) bases in the complement are characterized by 

the total entanglement of every qupit. This "standard distribution" is inescapable for two qupits 
. ^ (of any p), where only product and generalized Bell bases are admissible MUB types. This and the 

^ ! following results generalize previous results for qubits Q, and qutrits 16], drawing particularly 

upon Ref. [l3]. With three qupits there are three MUB types, and a number of combinations (p+2) 
are possible in full complements. With = 4, there are 6 MUB types for p = 2, but new MUB 
types become possible with larger p, and these are essential to realizing full complements. With 
this example, we argue that new MUB types, showing new entanglement characteristics, should 
enter with every step in N, and when is a prime plus 1, also at critical p values, p = N — 1. 
Such MUBs should play critical roles in filling complements. 

PACS numbers: 03.65.-w, 03.65.Aa, 03.65.Ta, 03.65.Ud. 



1 



I. INTRODUCTION 



raphy 



Mutual 



ly unbiased basis sets are known to provide an optimal basis for quantum tomog- 



2|, to play key roles in quantum cryptography 



solving the mean king problem in prime power dimensions 



-|6|, and to be instrumental in 
7| . The generalized Pauli opera- 



tors associated with MUB's include the stabilizers of quantum error correcting codes j8l-10|. 



and serve as entanglement witnesses 



for the MUB states. Of interest for the foundations 



of quantum physics, the MUB concept sharpens the concept of complementarity |l2l . 



and raises the question of existence in composite dimensions. An excellent comprehensive 
review of MUBs has recently appeared [l^ . 

We deal here with Hilbert spaces of prime power dimensions {d = p^), where d + 1 
MUBs are known to exist ^ . This is both the largest possible number, and also the number 
required for a complete operator basis (in representing the density matrix, for example). 
So, while each MUB is a complete orthonormal basis in the Hilbert space, the set of d + 1 
MUBs is a complete (nonorthogonal) basis in the space of all operators, which has dimension 
(p = p2Af Regarding terminology, to avoid reference to a "complete set of complete sets," 
and prompted by the fact that different MUBs (or the observable sets associated with them) 
are maximally complementary |l2|, I will use the term "full complement," or sometimes just 
"complement," to denote the set of all d + 1 MUBs. Partial MUB sets have been discussed 
in connection with composite dimensions and referred to as "constellations" 15 1. 

The natural systems to which MUBs apply consist of p-state objects (qupits). In such 
systems, while MUB complements exhibit only a single entanglement type for N = 2 (and 
all p), the number of distinct types proliferates with increasing A^. The variety is illustrated 
in a number of recent discussions, mostly on multiple qubit systems but also multiple qutrit 
systems g], [i3 
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2Q|. In particular, a systematic study by Romero and collaborators 
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illustrates a broad range of entanglement patterns that occur naturally in a construction 
scheme for full MUB complements. Such complements are catalogued for up to 4 qubits. 
Wiesniak and collaborators [l^ have developed a construction scheme aimed at experimental 
implementation and discussed the total entanglement content of full MUB complements of 
bipartite systems. 

With the general MUB problem in mind, our purpose here is to develop a general frame- 
work, independent of construction schemes, for exploring MUB entanglement patterns for 
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all p and N. The project begins by proving three general theorems (the "rules" ) that un- 
derlie and lead quickly to an array of more specific results. Many of the latter apply to 
all p, but are A'"-specific, as each step in introduces further richness. All results refer 
to one of two levels - that of individual MUBs and that of full complements. At the indi- 
vidual level, MUB types are characterized by first specifying the separation pattern - How 
many, and how big, are the irreducible subsets of qupits defined by the factorization of the 
wavefunction? - and next, by describing the entanglement pattern - What is the nature of 
the entanglement within each irreducible subset? At the level of the full complement, we 
ask about the possible MUB distributions - What combinations of MUB types can coexist 
within full complements. At the first level, we will show that all conceivable separation pat- 
terns are possible, and we will show with examples how to describe the entanglement within 
the nonseparable factors. At the level of the full complement, we will show how to deduce 
constraints on the possible MUB distributions. For N — 2 and 3, surprisingly, the general 
global constraints mentioned in the abstract suffice to determine all MUB distributions for 
all p. The N — 4 case is considerably more complex and requires the derivation of more 
detailed constraint equations. 

Let us begin with a review of basic concepts and notation in Section II. In Section III 
we prove the three general theorems. These rules are applied in Section IV to obtain the 
entanglement patterns of individual MUBs, and to deduce constraints on their possible 
distributions within full complements, taking the N — 2 - 4 cases in turn. In Section V we 
summarize results and comment on unresolved questions. 

II. BACKGROUND CONCEPTS AND DEFINITIONS 

In Hilbcrt spaces of dimension d, two orthonormal bases {K and Q) are mutually unbiased 
if any state \K,k) in basis K has uniform probably of being found in any state \Q,q) in 
basis Q; that is, if 



Thus, measurements in the two bases provide no redundant information. Since measure- 
ments in any basis provide d—1 independent probabilities, and since d^ — 1 real parameters 
are needed to determine an unknown quantum state (its density matrix p), it follows that 
d + 1 MUB's are required. In this way the MUB projectors form a complete nonorthog- 




(1) 
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onal basis in operator space. This required number of MUBs is (only) known to exist in 
power-of-prime dimensions. 



There is an intimate connection between MUBs and generahzed Pauh operators (herea: 
called simply "Pauli operators") which underlies several construction schemes (Ref. 



ter 



provides a comprehensive listing 2l|). These operators are conventionally written in the 
form of a tensor product, 

whose factors, acting on individual qupits, are powers of the generalized (p x -p) Pauli ma- 
trices, 

Z = ^\k)J'l,k\ and X = ^|A; + 1)(A;|, (3) 

fc=0 fc=0 

where tu = e^'^*/^, and X is the raising operator of Z . The powers n and m are p-nary 
numbers, eg^ n = {rii, ...niy), whose digits take the values 0,1,. ..,p — 1. Thus, there are p'^^ 
operators On,m (including the identity X = Oo,o); which make up a complete and orthonormal 
basis in operator space (with the trace operation as inner product). The desired connection 



with MUB's is described in Ref. j22|]: The On^m partition into p + 1 internally-commuting 
subsets, each consisting of p — 1 traceless operators (excluding X). The corresponding eigen- 
bases then form a complete complement of MUB's. 

The above are standard definitions and conventions. It will be useful to adopt a couple of 
more special conventions for use throughout this paper. First, the operator set On,m does not 
form a group, because multiplication generates irreducible phase factors. However, for odd p 
the set On,m ® (1? ...,uj^^^) does form a group, of order p'^^'^^, and for p = 2 the analogous 
set (9„,m® (±1, ±i) forms a grou p of order 2^^+^. These are called discrete Heisenberg-Weyl, 



or generalized Pauli groups |23|. We shall not make direct use of them, but we shall take 
advantage of the freedom to redefine the phases of the C„,m in the original set: We choose 
phases so that the compatible subsets form groups, and we call these compatibility groups. 
They are all isomorphic to those consisting of and Z"^, each of which is generated by 
the N independent elements, X = (Xi, ...,X^) and Z = {Zi, Z^), respectively. Thus, to 
construct another compatibility group, we may choose a generator set G = {Gi, ...,Gn) that 
consists of any N elements in the original compatible subset that do not form a subgroup. 
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and write the resulting group elements as 



G" = G;'^G^^..G7. (4) 

Thus, all of the compatibility groups are representations of the same group - the abelian 
group of order generated by elements. A simple example of a compatibility group so 
generated is 

= (XiZi)"i(X2Z2)"^..(XjvZ^)"-. (5) 

Note that phase factors are introduced with respect to the original Pauli operators because, 
eg, {X,Z,f = u'^^XfZl 

The generator set G, by itself, completely determines the states of the basis [G) in 
the Hilbert space, through the eigenvalue equations Gi\G,k) = u^^\G,k), where k = 
{ki, k2, kisi) is a p-nary representation of the state index k. The eigenvalues of a gen- 
eral group element are then given by 

=cj"-'=|G',A;), (6) 

where n ■ k = riiki + 712^2 + ■■• + rij^ik^, and the spectral representation of G" is therefore 
just the Fourier transform 2^ 

= ^ \G, k)Lo''-\G, k\ = J2 ^"^V{G, k), (7) 

k k 

where V{G, k) is the projector onto state k in basis G. This MUB projector is then given 
by the inverse transform, 

p(G,A;) =p-^^a;-"-^G". (8) 

n 

The existence of these simple transform relationships between every compatibility group and 
its corresponding MUB projector set is a consequence of defining the former to be a group. 
The only remaining arbitrary phases are those of the generators. 



III. GENERAL RESULTS ON ENTANGLEMENT 

In this section we establish the general rules that will form the basis for the rest of the 
work. For ease of reference and completeness I will state and prove these results as three 
separate numbered theorems. For transparency, here, in plain English, is what they will 
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say about MUB states: (I) A given qupit is perfectly pure or totally entangled, (II) The 
distribution of one-qupit operator factors in the compatibility group correlates with this 
purity,..., and (III) In any full MUB complement, every qupit appears pure p+1 times, and 
totally entangled — p times. 

These theorems and the results that follow from them rely on the assumption that MUB 
states are eigenstates of Pauli operators. While this is restrictive for individual MUB pairs, 
it is not restrictive for known MUB complements or known construction schemes 2l|, |25| . 
allowing for unitary equivalence. An example may help to illustrate. Consider the standard 
basis in 4D, and another related to it by the unitary transformation Un,m = {i)"'"^/2 (where 
n,m = 0,1,2,3), which is not an eigenbasis of the Pauli operators of Eq. [2J The two bases 
are MU, but a full complement cannot be completed containing both of them. However, 
full complements can be found containing either basis without the other: Starting with the 
well-known full complement containing the standard basis, one could apply Un^m to each of 
its bases to obtain another full complement. The latter are not eigenbases of the original 
Pauli operators, but clearly they are eigenbases of transformed Pauli operators, which may 
be thought of as corresponding to redefined parts (and redefined quantization axes). The 
results of this paper then apply with reference to these redefined parts. Regarding the 
existence of a MUB complement outside of this equivalence - I believe that this question 



also remains unresolved 



25| . We will return to these points in the conclusions. 



As a brief preliminary, one-qupit states within the A^-qupit system are defined by the 
reduced density matrices. 

Pi = Tr«p, (9) 

where Tr*^*^ denotes the partial trace over states of all but the i-th qupit. Perfect purity 
means that pi = pf is a projector, while total impurity means that pi = X/p. One can define 
the purity of the state pi as 

Pi = {pTrpl-l)/{p-l), (10) 
which takes its extremal values, 1 and 0, in the respective cases. 

Theorem I: If the system is in a pure eigenstate of Pauli operators (a generator set G), 
then any individual qupit must exist in a state of either perfect purity, or total impurity, 
the same for all eigenstates of G. 

Proof: The generators produce a compatibility group, and the A^-qupit density matrix 
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representing a pure eigenstate, p = V{G, k), may be expanded as in Eq. [HI Considering now 
the Pauli matrix factors that act on just the ith qupit, the generator set G must fall into 
one of two categories: Only one Pauli matrix, say Zi (and possibly powers of it), appears in 
the generator set, or more than one appear (including, say, Xi and Yi), that are not powers 
of one another. Consider the latter case, which is simpler: Let Gi and G2 be generators 
that contain the factors Xi and Yi. No operator of the form f/jX (where Ui is any one-qupit 
Pauli matrix) commutes with both Gi and G2, and all such operators are thereby excluded 
from the compatibility group. As a result, the only operator with a nonvanishing partial 
trace Tr^*) is the global identity X. Since X enters the summation with the coefficient 
, and Tr*^*) produces a factor of p^^~^\ the reduced density matrix for the ith qupit is 

pi = Xi, (11) 

indicating that the ith qupit is totally impure. 

Now turn to the other case: If only Zi (and possibly powers) appear in the generator set, 
then only Zi and its powers can appear in the compatibility group (again refering only to 
those factors that act on the ith qupit. Since the "one-body" operators Z"'X commute with 
all of these, they must belong to the compatibility group. These one-body operators are the 
only ones that survive the partial trace. Since each of them enters the summation (Eq. [8]) 
with coefficient and since Tr^*-* produces a factor of in each term, we find 

in this case that 

=p-i J^u;-"''^^^^ = \Z„h){Z„k,\. (12) 

n 

This shows that pi is a projector onto the eigenstate of Zi whose eigenvalue is u^', that is, 

p^i = Pi and ZiPi = id^'pi. (13) 

This proof is independent of the choice of the eigenstate k = {ki...k^) in the basis G, and 
so clearly the ith qubit is perfectly pure for all eigenstates in this basis . 

Here is a related more detailed theorem on the distribution of one-qupit matrices associ- 
ated with a single qupit. 

Theorem II: In any compatibility group of A^-qupit Pauli operators, the distribution of 
one-qupit factors acting on the ith qupit must be one of two types: (i) Only a single Pauli 
matrix and its powers occur, and each power occurs an equal number {p^~^) of times, or 
(ii) every Pauli matrix occurs, and each occurs an equal number {p^^^) of times. 
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Proof: Consider any set G of N generators of the compatibility group. This set must be 
one of the two types considered in the foregoing proof: Suppose first that only one Pauli 
matrix (say Zi), and possibly powers of Zi appear. Let Gi be a generator containing Zi as a 
factor, and let G2, G^,..., G^ he the rest. Gi by itself generates a cyclic subgroup containing 
all powers of Z^. Then, Gi and G2 by themselves generate a subgroup of order p"^ in which, 
by virtue of the rearrangement theorem, every power of Zi appears p times (no matter which 
power of Zi is present in G2)- One may repeat this argument, multiplying the order of the 
subgroup by p at each stage, until the full compatibility group is generated, with each power 
of Zi being produced p^~^ times. 

In the other case, let Gi and G2 be generators containing the Xi and Yi factors, respec- 
tively. These two generators, by themselves, generate a subgroup of order in which every 
Pauli matrix factor Ui appears once and only once. (To see this, note that Xi and Yi, by 
themselves, generate the one-qupit Pauli group [26|, but since Gi and G2 commute, the mul- 
tiplicity of phase factors is absent.) Now, by including a third generator, G3, one generates 
a subgroup of order p^ in which, by the rearrangement theorem, each Pauli matrix factor 
appears p times. Repeating the process through Gn, one generates the full compatibility 
group with each Pauli matrix factor appearing times. 

The second result is particularly striking in light of the fact that the nature of the entan- 
glement of the ith qupit may vary widely, in the sense that its entanglement may be shared 
with any number of other qupits in the system. Nevertheless, only two kinds of Pauli matrix 
distributions, with the correponding purities, are possible. 

We use both of the foregoing theorems to deduce the total entanglement content - as 
measured by the one-qupit purities - of a full complement of MUB's. This total content is 
constrained by the requirement that the two types of one-qupit Pauli matrix distributions 
be consisent with the set of all Pauli operators, which must appear in the full complement. 
Theorem III: Within any full complement of p^ + 1 MUB's, every qupit is perfectly pure 
in p + 1 basis sets, and totally entangled in the remaining — p. 

Proof: Consider the ith qupit. Recall that the total number of Pauli operators (excluding 
X) is p^^ — 1, and that these exactly accommodate the + 1 compatibility groups containing 
p^ — 1 traceless operators each. Each Pauli matrix factor Ui appears in p^A^-a pg^^n operators, 
except for Xj which appears in p2Af-2 _ ^ because we are not counting X in the individual 
groups. This number must equal the sum of Xj factors appearing in all of the compatibility 
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groups. According to the previous theorem, there are p^~^ — 1 such factors in compatibihty 
groups in which the ith qupit is pure, and p^^^ — 1 such factors in all other compatibility 
groups. If is the number of compatibility groups (or basis sets) in which it is pure, then, 
in order to account for all Xj factors, we must have 

p'""-' - 1 = i^sip""-' - 1) + (P^ + 1 - i^s)ip''-' - !)• (14) 
Solving this equation, we find the number of basis sets in which the ith qupit is pure, 

us = p + l, (15) 
and consequently, the number of basis sets in which it is totally entangled, 

^E = P^ -P- (16) 

The following corollary arises when all qupits take their pure states simultaneously: 
Corollary: The maximum number of product MUBs is p+ 1, and in any MUB complement 
where this number is realized, all of the remaining MUBs {p^ —p) must be totally entangled 
(in the sense that every qupit is totally entangled) [27|. This is the standard distribution. 

Note that the probability of finding the ith qupit pure in a MUB state picked at random 
from any full complement is equal to the averaged purity (Eq. [TOl) . 

J^s ^ P + 1 



(P-) = ' = ^ (17) 



which vanishes exponentially with A^. 



IV. ENTANGLEMENT PATTERNS AND THEIR STOICHIOMETRIES 

We discuss the N = 2 — A cases in turn. The first two are simpler, and we find that 
Theorems I and III are sufficient to determine all possible MUB distributions, although II 
provides useful insights. With = 4, we require Theorem II in deriving more detailed 
constraints that apply to individual qupits. 

bipartite systems 

Clearly, if one qupit is pure, then so must be the other. In light of Theorem I, then, 
both purities must be unity, or both zero. Because these purities coincide, the corollary of 
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Theorem III applies: There are p+1 product bases and —p totally entangled bases - the 
standard distribution is inevitable. 

We shall refer to all of the entangled bases as generalized Bell bases, because they share 
the common property that their compatibility groups consist solely of two-body operators, 



ie, those containing no factors [28[. To see the consequences of this, write one of the two 
generators as Gi = UV. The most general eigenstates of Gi may then be written as p-term 
expansions in the product basis of IV and UI, 

|^) = ^$^C^fc|A;)„|g-A;)„, (18) 

k 

where the eigenvalues of UV are u'^ and the coefficients Gk are determined by the other 
generator, call it G2 = ST. Commutativity demands that both S ^ U and T ^ V, so G2 
induces cyclic permutations (of order p) in the product states \k)y\q — k)^. Therefore the G^ 
are unimodular, and the p eigenvalues (w^) of G2 are nondegenerate, like those of Gi. This 
confirms explicitly what we know from Theorem I - namely, that measurements of one-qupit 
properties {eg, IV or UI) must produce random distributions over all possible outcomes. 
The generalized Bell states defined above are contained within a broader class definitions 



given elsewhere [29|, |30[ . The more restrictive definition given here - defining classes of states 
by the Pauli operators of which they form eigenbases - applies nonetheless to all MUBs 
that are compatible with known full complements, and we shall employ such definitions 
throughout this work as we proceed to larger N. 

We note for future reference that the precise form of the product state expansion f lTSj) 
depends on the choice of basis. A bad choice would require a p^-term expansion, but even a 
good choice could look slightly different. For example, if eigenstates of UV~'^ were expanded 
in the same product basis used in Eq. [181 one would find sums of + k)y. 

As a final note on Bell states, our working definition may be given in words alone: A 
generalized Bell state is any totally entangled two-qupit eigenstate of Pauli operators (since 
total entanglement requires that the two Pauli operators be of the form UV and ST). 

tripartite systems 

The standard MUB complement has p+1 product bases and p^ — p totally entangled 
bases. We shall refer to all of these totally entangled bases as generalized GHZ, or G-bases, 
because they have common properties describable as follows: 
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Let us first illustrate with a specific example that generalizes a standard choice of gener- 
ators for qubits |3l|, 



G = {Gi,G2,G3) = {XXY, XYX, YXX), (19) 

to arbitrary p. To identify an optimal product basis for an expansion, replace the latter two 
generators by G2 = G2Gi^ and G'^ = G^Gi^. Recalling the usual definition Yi = XiZi on 
the ith. qupit (modulo possible phase factors), the result is 

G' = {XXY, IZZ-\ ZIZ-^). (20) 

Clearly the most general joint eigenstates of G2 and G3 are p-term expansions in the standard 
basis, 

|V^) = ^J]C7fe|A; + g)|fc + r)|A;), (21) 

where and u'^ are the eigenvalues of and Gg, respectively, and the G^ are determined 
by Gi. The Gk are again unimodular because Gi generates a cyclic group of order p, of 
which the p product states form a basis. This again illustrates the randomness of one-qupit 
properties in totally entangled states. 

To demonstrate the commonality of all totally entangled three-qupit bases, we note that 
at least one generator must be a three-body operator (having no Ik factors), which we write 
in complete generality as Gi = UVW. Now, according to Theorem II, the inverse of each 
factor occurs p times in the compatibility group, once with the inverse of Gi itself, and 



p — 1 times in other three-body operators in which it is the only inverse (footnote 28|). 
Choosing two from the latter category, one containing and the other containing V~^, 
and multiplying Gi by each in turn, we obtain the generator set 

G = {UVW, IBG, AIG), (22) 

where compatibility requires that G is common to G2 and G3 as indicated. Clearly, A, B, 
and G define the product basis for the p-term expansions, 

1^) = ^ I] Ck\q - k)a\r - k),\k)„ (23) 

k 

and each of A, B, and G must differ from corresponding factors that appear in three-body 
operators of the compatibility group. In other words, every three-body operator in the 
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compatibility group induces cyclic permutations of the states composing the product basis. 
The similarity of generator sets shows that all totally entangled three-qupit bases have p- 
term expansions in some special product basis, and that all of their compatibility groups (of 
the same and p) have the same numbers of three-body and two-body operators. 

Again, a purely verbal definition is possible: A generalized GHZ state is any totally 
entangled 3-qupit eigenstate of Pauli operators. A general statement for A^ > 4 is possible 
but less categorical. 

The new aspect of MUBs that enters with A^ = 3 is the appearance of a third (nonstan- 
dard) MUB type, and with it, the possibility of composing a full complement with varying 
combinations. The third type is biseparable, and thereby nonsymmetric with respect to 
qupits - one qupit separates, leaving the other two in a Bell state. We shall refer to these as 
"separable-Bell" bases, with the shorthand notation SB (or SiB if we wish to identify the 



pure qupit). Such MUB bases are known for p = 2 and 3 (Refs. [l3|, ll6|), and to describe 
them for arbitrary p, we consider a generator set 

SiB = {I I A, UVI, ST I), (24) 

where UV and ST are commuting two-body operators acting on qupits 2 and 3. The 
joint eigenstates of this set may be written as 

\S,B:k,q,p) = \A,:k)\B2,^:q,p), (25) 

which describes qupit 1 in the fcth eigenstate of A, and qupits 2 and 3 in the Bell state 
denoted by the eigenvalues q and p of UV and ST, respectively. Similarly, the compatibility 
group of SiB is a tensor product of that associated with qupit 1 (Ji, Ai, A\~^) and that 
of the Bell basis of qupits 2 and 3. The tensor product is a common characteristic of all 
separable MUBs, and the eigenstates of a particular MUB all have the same character - 
the separation pattern involves the same entangled subsets of qupits, and the nature of the 
entanglement within each subset is the same. 

The three MUB types discussed above, including the three variations of the SB bases, 
exhaust all of the possibilities for three qupits. 

The remaining question now is, what combinations the three types of bases may appear 
in the full complement? One can answer this question simply by conserving the number of 
pure qupits while conserving the number of basis sets. We then find that we can remove a 
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TABLE I: Numbers of product, separable-Bell, and GHZ bases coexisting for three particles. 

single product basis (11) while adding three SB and removing two G bases: 

n + 2G^ 3SB (26) 

Table I shows the possibilities for three particles with any p. The cases of p = 2 and 
3 dramatize the role of totally entangled states with increasing dimension of the Hilbert 
space. In fact, case (a), dimension d — 8, is the only multiparticle MUB dimension in which 
a complement can be found with no totally entangled bases. And more typically, a majority 
of MUBs are totally entangled: In case (b) at least 4/7 of all bases are G bases, and even 
for two qutrits, 6 of the 10 bases are Bell bases. For = 3 and general p, the minimum 
number of G bases is given by Nmin{G) — p^ — 3p — 2, an ever-increasing fraction of the 
total number of bases as p increases. 

It is noteworthy that SB bases can be introduced only in steps of three, reflecting the 
condition that the three variations SiB must balance in the full complement, since the other 
MUB types arc symmetric with respect to permutations of qupits. This condition follows 
from the conservation of pure states for each qupit separately. 

quadrapetrtite systems 

The N — A case is more complex in a number of respects. Most importantly, new MUB 
types enter with increasing p. But even with p — 2, the number of distinct MUB types 
exceeds the number of separation patterns. Figure 1 shows the five separation patterns that 
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characterize all p, and lists seven MUB types, six of which account for all p = 2 options, and a 
seventh which represents, but is not exhaustive for p > 3. Let us first discuss the MUB types 
for general p, and later specialize to particular cases for constraints and stoichiometries. 



n 



SSB (six subtypes) 



SG(-^^ (four subtypes) 



BB (three subtypes) 



QC^) CC^) PC^) 



FIG. 1: Seven MUB types listed with 5 separation patterns for 4 qupits. 

The separable MUBs' compatibility groups are tensor products of those of their con- 
stituent MUBs, and their generator sets may be constructed accordingly. With SG^^\ a 
single generator is associated with the separating particle (for example II I A), while three 
generators (each of the form UVWI or their alternatives) are associated with the three 
particles forming G'-^^ states. There are four variations on this pattern corresponding to the 
choices of the separating particle. In the BB case, one could pick two generators of the form 
IIUV, and two of the form STII. There are three variations on this separation pattern, 
corresponding on the three ways of picking the two entangled pairs, as compared with six 
variations on the SSB pattern from the six ways of picking a single entangled pair. 

Let us discuss the nonseparable bases in somewhat more detail beginning with four- 
particle GHZ bases (G*^^-'). These are straightforward generalizations of the three-particle 
bases G^^\ and a standard generator set 3l| consists of the four operators 



G^^) = {XXXY, XXYX, XYXX, YXXX). 



(27) 



From an alternative generator set, {XXXY, IIZZ~^, IZIZ~^, ZIIZ~^), it is apparent 
that eigenstates may again be written as superpositions of p product states in the standard 
basis. A more general characterization of GHZ states is provided in the Appendix. 

Cluster bases {C^"^^) were introduced in connection with measurement-based, one-way 
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quantum computation 32|, and in fact both cluster and GHZ states are special cases of a 



broad class of A^-qubit states, called graph states, which form the basis of this [33|. Ref. 20 1 
has shown that graph states may be classified in terms of curves in phase space, which 
provides a further connection with the MUB problem. Cluster bases are defined here, for 
all p, by generator sets of which a standard example, introduced for the qubit case \22l, is 

= (XZXI, ZXIX, XIX Z, IX ZX). (28) 

Cluster states have stronger entanglement links between smaller groupings of particles, mak- 
ing their entanglement more .bust a,a„.t decoKe..ence Q than GHZ enta„,.e.e., wKich 

is shared equally among all particles. This is reflected in the fact that C^^'^ has only two 
2-body operators in its compatibility group, as compared with three in the G*-^-* case. For 
this reason, its generator set can only be simplified to {XZXI , ZIZ~^I ^ IZ^^IZ^ IX ZX), 
and as a result, the eigenstate expansions can be reduced to no less than terms in the 
standard basis. A general characterization of C*^^^ accompanies that of G*^^^ in the Appendix, 
which then goes on to show that these, together with the four separable bases, exhaust all 
MUB possibilities for four qubits. 

As a final example, I have found that a new type of basis, one that has no counterpart 
for qubits, is necessary for the existence of full MUB complements when p > 2, for reasons 
that will become apparent. A generator set giving rise to such a basis is 

= (ZXYW, XZWY, WYXZ, YWZX), (29) 

where standard definitions Y = XZ and W = XZ^ are followed. The essential point is that 
the generators are tensor products of four noncommuting one-body matrices, which rules out 
qubits, but makes possible the elimination of 2-body operators from the compatibility groups 
for p > 3. Less essential is that the four generators are related by pairwise permutations 
of operators (hence the notation P^'^^). The eigenstates have Bell correlations between all 
pairs of particles, not just the chosen pairs as in BB states. So, unlike cluster or BB states, 
the entanglement is shared equally among all four particles, but unlike GHZ states, the 
entanglement is robust. One can perform measurements on any two particles, in any two 
different bases, and produce a Bell state of the other two. 

To show that the P^^^ basis does not exhaust the possibilities for p > 3, we mention an- 
other generator set involving cyclic permutations, {ZXYW, XYW Z~^ , YW ZX, W Z~^XY), 
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where the factors are inserted for compatibihty. The corresponding basis could play a 
role similar to that of P'-^^ in filling MUB complements for p > 5, although it turns out to 
be relatively inconsequential when p = 3. In any case, since it would needlessly complicate 
the discussion without changing our conclusions, we exclude this example from the analysis. 

It is interesting to note in passing, that despite the differences in appearance among 
the generator sets of the four (five) totally entangled bases, the total numbers of factors 
appearing in their compatibility groups must be the same, namely 4 (p^ — 1), in accordance 
with Theorem II. This has consequences for stoichiometry, in particular for the standard 
distributions, and justifies classifying BB bases as totally entangled. 

Let us now turn to questions of stoichiometry. While in previous cases we were able 
to deduce the allowed entanglement patterns from global constraints alone (those involving 
total numbers of pure and entangled qupits in MUB complements), with > 4 this is no 
longer the case. The existence of multiple totally entangled basis types requires that we 
consider more microscopic constraints associated with the distributions of 1^ factors, as was 



done for qubits in Ref. [iTI]. To this end, we define a quantity that is capable of distinguishing 
among all MUB types under consideration. 

The "ra-body profile" of a particular MUB is the distribution of n-body operators (n = 1, 
2, A^) in its compatibility group, where (as implied earlier) n-body operators are those 
with N — n identity factors, Ik- This distribution is normalized to the total number of 
operators in the compatibility group, — 1. Examples of ra-body profiles are given in Table 
II, where we include the N = 2 and 3 cases both for comparison with = 4, and also 
to show how global information is recovered. The number of operators in each category, 
summed over all MUBs, must equal the numbers listed at the bottom of each column. The 
latter represent the n-hodj profile of the set of all Pauli operators, and are thus independent 
of the particular MUB choices. They are determined by generating all of the Pauli operators 
as expansions in the tensor products, 

(/i,Zi,Xi,...,XiZf-')) ® ... ® {In,Zn,Xn,...,XnZ^^-'^). (30) 

Thus, the total number of n-body operators is {'^){p — 1)", as shown. The condition that 
the MUB sums equal these bottom lines, column by column, provides (A^ — 1) independent 
constraint equations. (In the exceptional case of A^ = 2, both equations are independent.) 
It is immediately apparent from all of the first columns that the maximum number of 11 
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(b) = 3 

1-body 2-body 3-body 

n 3(p-l) 3(p-l)2 

SB p-1 

G 3(p-l) (p-l)2(p + 2) 

lOl 3(p2 - 1) 3(p2 - 1)2 (p2 _ 1)3 



(c) N = A 




1-body 


2-body 


3-body 




4-body 


n 


4{p - 1) 


6(p- 1)2 


4(p - 1)3 




(p-l)'i 




2(p-l) 


2p(p - 1) 


2(p-l)(p2-i) 


(P 


-l)3(p+l) 




P-1 


3(p-l) 


(p-l)2(p + 5) 


(P 


-l)'(p + 2) 


BB 





2(p2 - 1) 







(p2 - 1)2 







6(p-l) 


4(p-l)(p-2) 


4 

p - 


- 4p2 + 6p - 3 


C(4) 





2(p-l) 


4p(p — 1) 




- 4p2 + 2p + 1 


p(4) 








4(p2 - 1) 


p 


4 - 4p2 + 3 


all 


Ms? - 1) 


6(i?2 - 1)2 


4(p2 - 1)3 




(p' - 1)^ 



TABLE II: n-body profiles for all MUBs under discussion. The bottom lines ("all") are the n-body 
profiles of the set of all Pauli operators. 

bases is always given by p+ 1, the number that defines the standard MUB complement. Part 
(a) confirms that this is the only choice for N = 2, and its second column then determines 
the number of Bell bases (p2 — p), in accordance with the required total number of MUBs. 
Part (b) reproduces all = 3 results, as were summarized on Table I. The three columns 
provide three equations, but only two are linearly independent: The first column determines 
all possible combinations of IT and SB bases, and the second column then determines the 
number of G'-^) bases, which is again consistent with the required total number of MUBs, 
The third column provides no further constraint. 
Proceeding to the case of = 4, the calculation of the n-body profiles for the separable 
bases is straightforward, since their compatibility groups are tensor products of those whose 
profiles have already been calculated. The new nonseparable bases require more thought. 





(a) N = 


2 




1-body 


2-body 


n 


2(P-1) 


(p-l)2 


B 





p2-l 


all 


2(p2 - 1) 


(p2 - 1)2 
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1-body 2-body 3-body 4-body 



n 


4 


6 


4 


1 




2 


4 


6 


3 




1 


3 


7 


4 


BB 





6 





9 







6 





9 


C7(4) 





2 


8 


5 


all 


12 


54 


108 


81 



TABLE III: Specific n-body profile for four qubits. 

We found that the more symmetrical generator sets listed in Eqs. [2711291 were helpful in 
working out the profiles for general p. 

One can see by inspection of Table 11(c) that there is a qualitative difference between 
= 4 and the other cases. Consider just the first 6 MUB types, which represent all 
possibilities for p = 2. Looking at the 3-body factors in column (iii) , we can see that as 
p increases, the number of G^^^ and/or C^'^^ MUBs would have to increase as ~ in order 
to satisfy just Eq. (iii). But then they could not satisfy Eq. (ii), for they would produce 
too many two-body operators. Clearly, one eventually needs a basis which, like P^'^^ , has no 
two-body operators. This need makes itself felt already with p = 3, and becomes urgent with 
p = 5. With these differences in mind, let us consider the p = 2, 3 and 5 cases sequentially, 
to show how the general picture evolves with increasing p. 

four qubits 

The n-body profiles for p = 2 are shown on Table III. To explore stoichiometrics, con- 
sider the three equations (i, ii, and iii) represented by the first three columns, respectively. 
Equation (i), by itself, determines all possible combinations of the first three MUB types, 

4N{U) + 2N{S'^B) + N{SG^^^) = 12. (31) 

Next, notice that we can isolate the BB and G^^^ MUBs because of their simple profiles. 
Indeed, by simply adding (i) and (iii) we obtain the sum of all other MUBs, 

N{U) + N{S^B) + Ar(5G(3)) + NiC'^^^) = 15. (32) 
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{P = 


2) 


(P = 


3) 


(P = 


5) 


n 


3 





4 





6 










12 





16 





24 


BB 


2 


2 














C(4) 


12 


3 


72 


66 


360 


396 


p(4) 






6 





260 


206 


all 


17 


17 


82 


82 


626 


626 



TABLE IV: Examples of MUB distributions for four qupits with p = 2, 3, and 5. First, third, 
and fifth columns show standard distributions that maximize the numbers of totally entangled 
bases, while even columns show nonstandard distributions that minimize this number. Examples 
are chosen to minimize the number of P^^^ MUBs in all of the p = 3 and 5 cases. 

Since there are 17 MUBs in total we know immediately that 

N{BB) + N{G'^^'^)=2, (33) 

a result which also follows from 2(ii) + (iii) - (i), which reproduces the total number. 

There are 16 ways to satisfy Eq. [3ll with N{C'^^'') determined in each case by Eq. [321 For 
each of these combinations, there are 3 ways to satisfy Eq. [SSI for a total of 48 possible MUB 
distributions. To illustrate the range, a standard and a nonstandard distribution are shown 
in the two leftmost columns of Table IV. These examples are chosen to show the maximum 
and minimum numbers of the new (nonseparable) C^^'^ MUBs, which make up a majority 
of MUBs in 30 of the 48 possible distributions. The dominance of C^^^ MUBs is related to 
the large number of 3-body operators in their profile. Similar complements were found in 
Ref. |l7| through an explicit construction, except that the G^^-* MUBs were not produced, 
so that 16 combinations were obtained with 2 BB MUBs present in all of them. 

four qutrits 

The n-body profiles for the p = 3 case are shown in Table V. Again, the first column 
restricts the combinations of the first three MUB types, 

4A^(n) + 2N{S'^B) + N{SG^^'^) = 16. (34) 
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(a) p = 


3 






(b) p = 


5 






1-body 


2-body 


3-body 




1-body 


2-body 


3-body 


n 


8 


24 


32 


n 


16 


96 


256 


S^B 


4 


12 


32 


S^B 


8 


40 


192 




2 


6 


32 




4 


12 


160 


BB 





16 





BB 





48 










12 


8 


G(^) 





24 


48 


C(4) 





4 


24 







8 


80 


p(4) 








32 


p(4) 








96 


all 


32 


384 


2048 


all 


96 


3456 


55296 



TABLE V: The n-body profiles for four qutrits and four ququints. 

The first and second columns together [(ii)— 3(i)] restrict other combinations, 

4N{BB) + 3iV(G'(^)) + iV(C(^)) = 72, (35) 

and the inchision of the third column [(iii)+2(ii)— 6(i)] yields the total MUB count, 

Ar(n) + ... + Ar(p(4)) = 82. (36) 

There are 25 combinations of the first three MUB types that satisfy Eq. [311 as compared 
with 16 such combinations in the qubit case. But, in the absence of P^^^ MUBs, one cannot 
solve both Eqs. [35] and [36] for all of these combinations, and we find a total of only 11 MUB 
distributions. To trace the reasons, we subtract Eq. [35] from [36] and solve for P^^'^: 

iV(p(^)) = 10 + 3N{BB) + 2N{G^^^) - [N{U) + N{SSB) + N{SG^'^^)]. (37) 

Without P^^^ MUBs the left side vanishes, and there can be solutions only if the quantity in 
square brackets is 10 or larger. This condition fails for the standard distribution, for which 
this number is A^(n) = 4. In this case, the minimum number of P^^'^ MUBs is 6, as shown 
on Table IV. The other entry maximizes the quantity in square brackets at iV(5G'(3)) = 16. 
In both entries we then minimize the number of P'-^^ MUBs by maximizing the number of 
MUBs. One can increase the number of P^"^^ MUBs over these minima by adding BB 
and/or G^'^^ and subtracting C^^^ MUBs. Although its numbers can be small, the P(^) MUBs 
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play a critical role in maintaining the balance of 3-body operators (Table V) at no cost in 
two-body operators. 

With P*^^) MUBs included, the multiplicity of each of the 25 solutions of Eq. |21] is large 
(we estimate more than 200), for a total of probably more than 5000 solutions. We cannot 
argue that all of these solutions represent realizable MUB distributions, because we cannot 
rule out the possibility of more subtle constraints. Such concerns are beyond the scope of 
the present paper. 

four ququints 

Again consulting Table V for the p = 5 case, it is striking to see how three simple equations 
can again emerge from appropriate combinations. The first column gives us directly 

4N{U) + 2N{S^B) + NiSG^^')) = 24, (38) 

the combination [(iii)+8(ii)— 64(i)] relates the other four quantities, 

8N{BB) + 5N{G^^^) + 3Ar(C(^)) + 2N{P'^^^) = 1600, (39) 

and still another combination [(iii)+2(ii)— 22(i)] yields the total MUB count, 

A^(n) + ... + N{P^^^) = 626. (40) 

There are 49 combinations of the first three MUB types that satisfy Eq. [38l but in the 
absence of the P(^) MUBs, none of these admits solutions of Eqs. [39] and HQl To see how 
this situation arises, solve the latter two equations for N{P^^^) while eliminating N(C^^^): 

Ar(p(^)) = 278 + 5N{BB) + 2N{G^^^) - 3[N{U) + NiSG'^^^) + N{SSB)]. (41) 

The quantity in square brackets has minimum and maximum values of 6 (the standard 
distribution) and 24, as shown on Table V, corresponding to lower bounds on N(P'^'^^) of 
260 and 206, respectively. The latter is the absolute minimum number of P(^) MUBs m any 
full complement. Again, one can add P(^) MUBs by removing G^'^^ and adding BB and/or 
G(^) MUBs, so that pW can be the majority MUB type in some complements. While P^^^ 
is critical for both p = 3 and 5 cases, it plays a considerably more dominant role here. The 
underlying reason is that the ratio of the numbers of 3-body to 2-body operators increases 
considerably in going from p = 3 to 5, as shown in Table V. 
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We estimate the total number of solutions of Eqs. l38M0l to be in excess of 10^, but again, 
we cannot argue that all such solutions represent realizable MUB distributions, or provide 
a revised estimate, without a further study of possible constraints. 

The examples of this section have shown us that with every step in A^, and with some 
steps in p, full complements require not only those MUB types generated from smaller sys- 
tems, but also new, nonseparable MUB types that exhibit new entanglement characteristics 
inaccessible to smaller systems. In the step from = 2 to 3, G^^^ MUBs are required for 
the standard distribution, although a nonstandard distribution {SB only) is possible with 
p = 2. With the step to p = 3, no MUB complement exists without G^^\ In the step to 

= 4, the ^(4) MUBs are indespensible to all MUB distributions with p = 2. With the 
step to p = 3, the new P^^^ MUBs become possible, and they in turn make possible the 
standard distribution. At p = 5, the P^'^^ MUBs become indispensable to all distributions. 

Projecting to larger systems, the distinguishing feature of the P^^^ generator set is that 
a different (noncommuting) Pauli matrix factor is associated with each qupit. The number 
of such factors in general is p + 1, and when this is equal to the number of qupits, a new 
type of entanglement becomes possible. Thus we predict that when A^ is equal to any prime 
plus 1, then that prime {p^ = A^ — 1) is a critical value for the emergence of new entangled 
states as p is increased at fixed N. These states should play critical roles in filling MUB 
complements for p equal to or slightly greater than p^r. 

V. CONCLUSIONS AND OPEN QUESTIONS 

We have exploited the connections between MUBs and Pauli operators to develop a 
general framework for investigating both the entanglement properties of individual MUBs, 
and the combinations of such MUBs that can be found in full complements. We began by 
proving general theorems regarding MUBs as eigenbases of Pauli operators: We showed that 
the purities of individual qupits in such eigenbases must be either or 1, that the purity alone 
dictates the distribution of Pauli matrix factors (including J^) in the compatibility groups 
of these MUBs, and that every qupit must adopt these special purities the same number 
of times within any MUB complement: (p + 1) times pure, and (p^ — p) times totally 
entangled. An immediate corollary is that one may have at most p + 1 product MUBs in a 
full complement, and when one does, all remaining MUBs must be totally entangled. This 
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defines the standard distribution. 

Armed with these theorems and the general properties of Pauh operators, one quickly 
obtains more specific results: When N = 2, only product and generalized Bell bases are 
possible, for any p, and the standard MUB distribution is inevitable. With = 3, the 
unique totally entangled bases are generalized GHZ bases, but a third MUB type becomes 
possible, namely separable-Bell bases. This makes possible p+2 distinct MUB distributions. 
With N — A and p — 2 there are six MUB types, including two nonseparable bases and a 
third (BB) that is separable but totally entangled. There are 48 possible MUB distributions, 
with cluster bases making up the majority of MUBs in most of these. With N — A and 
larger p, further MUB types exist, and at least one such MUB type (P^^^) is essential to 
forming a standard MUB complement with p = 3, and to forming any MUB complement 
with p — 5. 

Several results have emerged in the course of working the above examples, and it seems 
useful to synthesize these in one place: (1) A MUB can exist in any separation pattern. 
(2) All states in a particular MUB have common separation and entanglement patterns - 
the generator set contains all information about the nature of the entanglement, while the 
eigenvalues specify the states. (3) Compatibility groups of separable bases are tensor prod- 
ucts of those of the nonseparable constituent bases. It follows that (4) within nonseparable 
groupings of qupits, those with two qupits must be in generalized Bell states, those with 
three qupits - generahzed GHZ states. Those with 4 qupits have the same broader array of 
options available to 4-qupit systems. 

Perhaps the most important lesson to be drawn from the present examples is that, al- 
though it is easy to construct MUBs from those found at lower A^, either as tensor products, 
or as larger- A^ counterparts such as G^^\ the more interesting challenge is to find the new 
nonseparable MUB types, with no counterparts at smaller A^ (or sometimes p), that make 
full complements possible. It may be a general feature that such MUBs tend to dominate 
MUB distributions near the A" and p values where they first emerge, only to be superceeded 
by other MUB types as the system size increases. In this sense, every A" is critical, but not 
every p. We predict that when A^ is a prime plus 1 (eg, N = 6, 8,...), there will be a critical 
value, Pat = A^ — 1, for the introduction of new entangled states that will play critical roles 
in MUB distributions. 

Glosing thoughts on the existence question: The intimate connection between MUBs 
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and entanglement for > 2 highlights the way in which all known MUB complements 
;ake advantage of the symmetry associated with equivalent parts, unique to dimensions 



34|- Theorem I, which makes no reference to dimension, cannot hold in (at least some) 
composite dimensions: In the simplest counterexample, 6 dimensions, the qubit can be 
totally entangled, but the qutrit cannot be. Yet, the import of Theorem I is that for all 
known MUB complements, there exists a factorization into parts (represented by some set 
of generalized Pauli operators), in terms of which Theorem I (and the others) hold. Thus, 
in addition to the existence question in composite dimensions, there is also an existence 
question in dimensions - Do MUB complements exist that violate Theorem I? Perhaps 
entanglement considerations such as the present ones will help in answering these persistent 
questions. 
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APPENDIX 

First, we write general definitions of G^^"^ and C^^^ bases in terms of generator sets, and 
then we show that these are the only possible nonseparable bases for four qubits. 

To generalization Eqs. [27] and [28] with maximum transparency, we follow the alternative 
forms written in the text and define G^^-* bases by the generator set 

G(^) = {ABII, AICI, AND, STUV), (42) 

where every four-body factor must differ from its two-body counterpart {S ^ A, etc.). 
The two-body operators provide the special product basis for the p-teim expansion. The 
generalization to A^-qupit GHZ states is apparent. 

The C^^^ bases are best defined in a similar way, although a bit less transparently because 
there are only two independent two-body operators, 

C(^) = {AICI, IBID, SBUI, ITCV). (43) 

Again, the two-body operators provide a product basis for the expansion, which in this 
case requires terms. Individual factors in the three-body operators must differ from 
corresponding factors in the two-body operators, except where their equality is explicit. 
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It should be noted that there are two variations on the C*^^^ generator set, corre- 
sponding to the other ways of pairing the two-body factors. One such variation is 
[ABII, IICD, STCI, IBUV). The three alternatives are mathematically equivalent, 
although one can make a physical distinction based on entanglement links between pairs. 
The stronger entanglement links in a system represented by Eq. |33] are between neighbors 
in the sequence (1-2-3-4-1). In the variation given above, the sequence is (1-3-2-4-1), and in 
the other possible variation it is (1-2-4-3-1). The various possibilities are not unphysical, as 
one can imagine unlike particles with tetrahedral coordination. 

Completeness for qubits 

We now argue that the two bases defined above are the only nonseparable options for 
qubits. We first argue that both are the unique nonseparable representatives of their respec- 
tive n-hodj profiles as shown on Table III. We then show that other profiles cannot exist 
for qubits. 

It is straightforward to verify that the G*-^^ generators produce six two-body operators 
involving the same factors, A - D. These exhaust all 12 of the available Ik factors (Theorem 
II demands three factors per qubit), so that all remaining operators are 4-body operators, 
as shown in the table. The only other basis that can share this profile is BB, in which the 
two-body operators have factors that do not commute individually. But, by virtue of this 
fact, the BB basis has four independent two-body operators, so that these can compose the 
generator set. The separability of the basis is then obvious. 

Turning to the C^^^ case, it is straightforward to show that the generators of Eq. |13] 
produce no further two-body operators beyond the two shown, so that remaining 1^ factors 
must appear with the 8 three-body operators. One might wonder whether a different basis 
could be found with the same profile by using four-body generators in place of the three- 
body generators. The answer is no - It is easy show that this would generate only G*-^-* or 
BB bases, depending upon whether one of the four-body generators shares factors with one 
of the two-body generators. 

The remaining point is to rule out other four-qubit profiles. It suffices to consider just 
the two-body operators, whose maximum number is six. We will show that the numbers 
4 and are impossible for qubits. The former case is very simple - it is impossible to 
find four commuting two-body operators that do not generate two more (and these will 
immediately identify themselves as belonging to either a G^^^ or BB compatibility group). 
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As to the latter case, assume that there are no two-body operators. Then all 12 of the 
Ik factors must appear in one-body operators, making three appearances on each qubit 
(Theorem II). Consider any two of these operators that have their 1^ factors on the same 
qubit. Commutativity demands that they have exactly one other factor in common, so that 
their product is a two-body operator. This forms a contradiction and shows that there is 
no profile without two-body operators. 

To briefly summarize the results of this appendix, the 6 MUB types hsted on Table 
III exhaust the possibihties for four qubits. The 5 corresponding A^-body profiles are also 
exhaustive; in particular, a P^^^-like profile does not exist for qubits. 
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